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The onset of equilibration of partons in nuclear collisions may lead to related trends 
in the centrality dependence of {p t ) ) p t fluctuations, and net charge fluctuations [Q], We 
extend the transport description of ref. [ lj to include radial flow. 

1. INTRODUCTION 

Liquids exhibit shorter correlation lengths than gases. Correlation and fluctuation 
measurements can therefore reveal the state - liquid or gaseous - of the matter produced in 
heavy ion collisions, but we must first sort out the ion-collision effects that also contribute. 
Experiments show that dynamic fluctuations of the transverse momentum are roughly 
independent of beam energy, but increase as centrality increases [ El HI IB] Local 
equilibration can explain this behavior [Q], but flow and jet explanations have also been 
suggested [ElEj- Certainly all effects play a role at some level. It is therefore important 
to describe all of them in a common framework. We take a modest step in this direction 
by extending the formulation in [ lj to account for contributions from radial flow. 

Dynamic fluctuations are obtained from the measured fluctuations by subtracting the 
statistical value expected in global equilibrium [ 2j. For particles of momenta pi and p 2 , 
dynamic multiplicity fluctuations are characterized by 



R = (N ) ~ ( N r - (N) 



^2 ^2 / rf Pl rf P2KPl,P2), (1) 

where (• • •) is the event average. This quantity depends only on the two-body correlation 
function r(pi, p 2 ) = iV(pi, P2)— iV(pi)iV(p2). It is obtained from the multiplicity variance 
by subtracting its Poisson value (N). For dynamic p t fluctuations one similarly finds 

where 5p ti = p ti — (p t ). CERES and STAR measure this observable, while PHENIX 
measures F pt « N(Sp t iSp t 2) /2cr 2 , where a 2 = (p 2 ) — (pt) 2 - 

Figure 1 shows the ratio (SptiSp^) 1 ^ 2 / (pt) for central collisions as a function of beam 
energy as compiled in [ 2 J from CERES and STAR data [El Ej- The data suggest that 
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Figure 1. Energy dependence of p t fluctua- 
tions from [ 2j. CERES and STAR data are 
from [ El E] . The solid line follows from (JBJ) . 



Figure 2. Centrality dependence of p t fluc- 
tuations [ 4j. Centrality is determined by 
the number of participants relative to the 
maximum value. 



the source of fluctuations is independent of energy, although uncertainties are large. This 
energy independence suggests that jets are not the leading contribution to p t fluctuations. 

2. THERMALIZATION 

In [ lj, the Boltzmann transport equation is used to show that thermalization alters 
the average transverse momentum according to 

(pt) = (pt}oS+(p t }e(l-S), (3) 

where S = is the probability that a particle escapes the collision volume without 
scattering. The initial value (pt)o is determined by the particle production mechanism. 
If the number of collisions Af is small, S ~ 1 — Af implies a random- walk-like increase of 
(p t ) relative to (pt)o- As centrality increases, the system lifetime and Af both increase, 
eventually to a point where local equilibrium is reached. Correspondingly, the survival 
probability S in Q decreases with increasing centrality. The average (p t ) increases for 
more central collisions to the point where local equilibrium is established. The behavior 
in events beyond that centrality depends on how the subsequent hydrodynamic evolution 
changes (p t ) e as the lifetime increases, as discussed in the next section. 

Dynamic fluctuations depend on two-body correlations and, correspondingly, are quadratic 
in S. In [ lj, Langevin noise is added to the Boltzmann equation to describe the fluctua- 
tions of the phase space distribution. If the initial correlations are not far from the local 
equilibrium value, we find 

{SptiSpta) = {SptiS Pt 2)oS 2 + (S Ptl S Pt2 ) e (l - S 2 ). (4) 

As before, the initial quantity (Sp t iSp t 2) is determined by the particle production mech- 
anism, while (SptiSp t 2)e describes the system near local equilibrium. 

Near local equilibrium, spatial correlations occur because the fluid is inhomogeneous - 
it is more likely to find particles where the density is high. The mean p t at each point is 
proportional to the temperature T(x), so that 

(5pn5p t2 )e~ J r(xi J x 2 )5r(xi)5r(x 2 ) J (5) 
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where r(xi, x 2 ) is the spatial correlation function and ST = T — (T). Our comparison to 
the STAR data computed in [ lj is shown in fig. 2. 



where R is the scaled variance (PJ. The dimensionless quantity F depends on the ratio 
of the correlation length to the transverse size. In [ lj, we took F and R to be energy 
independent. Therefore, if matter is locally equilibrated in the most central collisions, we 
would obtain the solid line in fig. 1. 

We do not expect F to vary strongly with energy unless the correlation length changes, 
e.g., from a liquid to a gas value. In addition, we expect the scaled variance R to vary 
primarily with the number of participants and, hence, to be roughly independent of energy. 
For comparison, the dashed curve in fig. 1 shows an ad hoc TV -1 / 2 dependence on the 
measured multiplicity. 



Thermalization describes the p t fluctuations measured by STAR [Hj in peripheral colli- 
sions as shown in fig. 2, but model calculations fall short of the data for the most central 
collisions. One explanation of this discrepancy is the appearance of radial collective flow 
in the nearly-thermalized central collisions. Radial flow was omitted in [ lj, which only 
addressed longitudinal expansion. Flow arises when pressure gradients within the fluid 
cause it to move outward with an average local velocity. This collective motion requires 
the fluid to be sufficiently near local equilibrium that each fluid cell can respond to the 
different pressures of nearby cells. 

To see how flow fits into our transport framework, consider a test particle traveling 
through the medium. As the system begins to thermalize, the test particle undergoes 
a random walk so that its root-mean-square p t increases. At the same time, all of the 
other particles are experiencing their own random walks in both position and momentum 
space. On average, the walk in position space directs particles radially outward with an 
average drift velocity taking along the test particle. Collective flow builds as more of 
the particles encountered by the test particle are moving radially outward. We expect 
this effect to be most important in central collisions, since the system lives longer and 
particles have more time to interact. 

A simple way to describe flow is to introduce an average radial parton drift velocity v. 
This drift causes a radial blue-shift for massless partons, so that the parton Boltzmann 
distribution is /(p) = exp(-£/T eff ), where T eff = T[(l + v)/(l - v)] 1 ' 2 . The average p t 
and its fluctuations near equilibrium then satisfy 




3. FLOW 




and 




(7) 



Observe that (SpuSp^) 1 ^ 2 / (pt) is independent of the flow velocity v\ fig. 1 is unaffected. 
We obtain the mean velocity as a function of centrality by fitting the (p t ) measured by 
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Figure 3. Centrality dependence of p t fluctu- 
ations with radial flow included. 



Figure 4. PHENIX data for central 
Au+Au for p t < p^max [ 5j. The solid 
curve includes flow, while dashed curve 
does not. 



PHENIX for pions while using the same (pt)° e as before [Ej. The net effect of including 
velocity in |Zj) is to compensate the longitudinal cooling. 

By including flow we describe the data for the full centrality range, see fig. 3. Flow 
effects are further exhibited by measuring fluctuations in the acceptance region 0.2 GeV < 
Pt < Pt, max as a function of p^ max following PHENIX [ 5j. Observe that © implies 
F pt « N(5p t i5p t2 ) /2a 2 oc NR(p t ) 2 /(j 2 ~ NR. The multiplicity N grows strongly as f^ max 
is increased, since N is the integral of exp(— Pt/T e ^) over the acceptance. To see that R is 
independent of p^ max , let e be the probability that a particle falls in the acceptance region 
and 1 — e be the probability that it is missed. The average number of detected particles 
is (N) e = e(N). For a binomial distribution, one has (N 2 ) e = e 2 (N 2 ) + e(l - e)(N). 
Equation (P) then implies that R(e) = i?, the value for full acceptance. 

The dashed curve in fig.0|is computed from (JOJ) at a freeze out temperature fixed by the 
measured (p t ). The solid curve includes flow through T e ^ and fits the data much better. 

In summary, we find that flow is an important aspect of the multiple scattering contri- 
bution to p t fluctuations in central collisions. The flow and thermalization contributions 
are essentially independent of beam energy in accord with data [Ej. Any jet contribution 
to this phenomena would be tightly constrained by these data [Ej. 
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